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ABSTRACT 

The fundamental parameters describing the coronal spectrum of an accreting black hole are the slope 
P of the power-law continuum and the energy at which it rolls over. Remarkably, this latter 
parameter can be accurately measured for values as high as 1 MeV by modeling the spectrum of X- 
rays reflected from a black hole accretion disk at energies below 100 keV. This is possible because the 
details in the reflection spectrum, rich in fluorescent lines and other atomic features, are very sensitive 
to the spectral shape of the hardest coronal radiation illuminating the disk. We show that by htting 
simultaneous NuSTAR (3-79 keV) and low-energy (e.g., Suzaku) data with the most recent version of 
our reflection model relxill one can obtain reasonable constraints on E^ut at energies from tens of 
keV up to 1 MeV, for a source as faint as 1 mCrab in a 100 ks observation. 
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1. INTRODUCTION 

X-ray reflection spectroscopy is currently our most ef¬ 
fective tool for probing accretion processes and space- 
time near the supermassive black holes that power ac¬ 
tive galactic nuclei (AGN). In the current paradigm, the 


power-law continuum is produced either in a central hot 

corona (e.g., Shakura & Sunvaev||1973 

Haardt||1993D or 

at the base of a jet (e.g., Alatt et al.| 

TF9 

2 iVlarkotf et af. 

2005). A fraction of this radiation ifl 
tion disk and is reprocessed into a 

uminates the accre- 
rich reflection spec- 


trum consisting of fluorescent lines and other features, 
which provides detailed information on the composition 
and ionization state of the accretion disk and constraints 
on the structure of the corona. Importantly, reflection 
spectroscopy of AGN is the way to infer the spin of their 
black holes via modeling key spectral features, such as 
the Fe K line (~ 6.4-6.9 keV), the Fe K edge (~ 7- 
8 keV), and the Compton hump 20-40 keV). Obser¬ 
vational data suggest that in many cases the illuminating 
power-law continuum cuts off abr uptly in the hard X-ray 
band (e.g., Zdziarski et al.|[2000|. Most X-ray detectors 


are unable to accurately constrain this cutoff energy be¬ 
cause they either have limited sensitivity at high-energies 
or none at all. Thus, reflection models have been tradi¬ 
tionally simplified by adopting a fixed c utoff energy at a 


reasonable value, e.g., Er.,,, = 300 keV (Ross & Fabian 

2005 Garcia & Kallman||2010|. _ 

'i'he 2012 launch of NuS'l'AR (Harrison et al. 20131 
has revolutionized reflection spectroscopy because its in¬ 
struments provide both very low background and superb 
sensitivity over the band 3-7 9 keV, which captures all the 
key reflection features (e.g., Risaliti et al. 2013 Walton 


^ Harvard-Smithsonian Center for 


Astro^hysics^ 60 ^_Gard^ 

lerfSheaa.eta. harvard.edu, 


""ieTTiHcla.harvard.edu. isteTrier^TieaaTTaTiarvarfrfinTT 


^ Dr. Karl Hemeis-Observatory and Erlangen Centre for As- 

troparticie Fhysics, !Sternwartstr 7, 95049 B amberg, Germai^ 


mas'nauserHHernwarte'mi^rTangen.de" 


^ Institute for Astrophysical Kesearch, Boston University. 

^25 CJommonwealth Avenue, Boston, MA 02215, UsA; 


keckmObu.edu 


et al.||2013 jParker et al.||2014|. The quality of the NuS- 
7AA data demand reflection models that properly treat 
the high-energy cutoff as a fit parameter. We have de¬ 
veloped advanced relativistic models of X-ray reflection 
from ioni zed accretion disks; the most recent version is 
relxill (Garcia et al. 2014). These models include a 


rich atomic database, tuff treatment of the angular dis¬ 
tribution of the reflected radiation, and an improved ge¬ 
ometrical model of the illuminating coronal source. Fur¬ 
thermore, the most recent version of relxill also in¬ 
cludes the high-energy cutoff Ecut as a ht parameter. 
Likewise, E^-at is included as a free parameter in t he most 
recent version of the reflection code reflionx (Ross & 
Fabian|[200^ . 


in the standard picture, the power law spectrum is gen¬ 
erated in the corona by Compton up-scattering of ther¬ 
mal disk photons. The shape of the spectrum depends 
on the geom etry, temperature, and opti cal depth of the 
corona (e.g., Rybicki & Lightman 1979). In particular, 
the fractional energy change of a scattered photon (in 
the non-relativistic limit) 


AA 

IT 


dfcTe 

TTleC^ 


E 


( 1 ) 


implies a cutoff of the power-law near 4fcTe because 
the energy gained by the photon cannot e xceed the en¬ 
ergy of the electron (e.g., Titarchuk 1994). In practice, 
Acut~ 2 — 3fcTe due to dispersion, geome trical, and rela¬ 


tivistic effects (e.g., Petrucci et al. 


2001). Thus, an em¬ 


pirical estimate of L’cut provides direct information on 
the temperature Tg of the coronal electrons. 

The recent data obtained with NuSTAR has spot- 


lighted Ecut (e.g.. 

Brenneman et al. 

2014 

Ballantyne 

et al.|2014| 

Marinucci et ai.|2014 

Halo 

kovic et ai.|2015 ). 


these studies, the coverage of NuSTAR to 
79 keV was sufficient to directly observe the roll over 
of the power-law component {Ecut ~ 50 — 200 keV). 
However, fitting a NuSTAR s pectrum of I MGC 5506 with 
relativistic reflection models, |Matt et al. (2015) surpris- 

in fitting a 


ingly found Ecut=720ljgg keV. Likewise, 
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Suzaku+NuSTAR spectrum of NGC 4151 [Keck et al7| 
(2015) found that a value of i?cut near 1 MeV was sta- 


tistically required. 

In this Letter, we demonstrate the remarkable ability 
of reflection models to provide strong constraints on E^ut 
up to several hundred keV, values that greatly exceed the 
limit of the detector bandpass (e.g., 79 keV for NuSTAR). 
While not detected directly, this high-energy portion of 
the spectrum conditions the ionization state and struc¬ 
ture of the disk atmosphere, which in turn modifies the 
observed reflection features. Thus, by properly modeling 
the reflection component at moderate energies, one can 
accurately estimate Lieut at extreme energies and thereby 
constrain the physical properties of the corona. 

2. MODELING THE REELECTION SPECTRUM 
The reflection model relxill is a fusion of our reflec- 


tion code xillver (Garcia &: KallmanI 


ian|20101 
20IU1 12(]' 


Garcia et al. 


2013) and relline (Dauser et al. 21)11) '5oT3), which is 
a general-relativistic ray tracing code. We hrst describe 
the atomic physics part xillver, which assumes plane- 
parallel geometry and idealizes the disk as a slab with a 
total optical depth of r = 10 and a constant density of 
n = 10^® cm“^. The illuminating source is assumed to 
have a power-law spectrum with a photon index T and 
an exponential cutoff at high energies, £lcut0 The in¬ 
tensity of the illumination is controlled by specifying the 
ionization parameter: 




47rF. 


( 2 ) 


where L"x is the net ionizing flux in the energy band 1- 
1000 Ry. A comprehensive grid of models has been calcu¬ 
lated and provided to the community for fitting observa¬ 
tional data that cover a wide range of model parameters. 
The most recent version of xillver includes models for 
several different values of the high-energy cutoff, namely, 
= 20,50,80,100,150,200,250,300,500,1000 keV. 
In Figure we show examples of models for these val¬ 
ues of Ecnt for two ionization parameters, log ^ = 1 and 
log C = 3. 

As n iustrated in Figure [T] and discussed in [Garcia et al.| 
(2013), Acut controls the intensity and shape of the il- 


lummating spectrum and thereby affects the global re¬ 
flection spectrum. The Compton hump is most dra¬ 
matically affected, its peak shifting to lower energies as 
Ecut decreases. Weaker effects are seen at lower energies 
(< 10 keV), which are mostly driven by changes in the 
ionization state of the gas. As Equation [^suggests, two 
models with the same ionization parameter (and hence 
the same E^) but different values of Acut will give rise 
to distinctive reflection spectra. This difference occurs 
because the illuminating spectrum with lower E^ut will 
have a significantly larger proportion of its photons con¬ 
centrated at lower energies, which favors photoelectric 
interactions over electron scattering relative to the spec¬ 
trum with high Ecut- 

Changes in E^^t and the consequent shift in the balance 
of photons in the hard and soft bands affects the ioniza¬ 
tion structure of the slab, as illustrated in the two right 

^ For simplicity, we use a phenomenological e-folded power-law 
model rat her than a proper mod el of thermal Comptonization such 
as that of |Zdzia rski^e|j Ln(|2003|. a rnodel whose accuracy has been 
demonstrated (e.g., jtiieninski et al.|1997| jWardzihski i 


et al.|2002|. 


panels of Figure On the one hand, for low ionization 
the lower values of Ecut (strong soft flux) produce hotter 
atmospheres; as a consequence, in this regime the tem¬ 
perature is mainly regulated by photoionization heating 
and recombination cooling. On the other hand, for high 
ionization the hotter atmospheres correspond to higher 
values of Ecut (strong hard flux); in this regime, the 
temperature is controlled by electron scattering of the 
abundant hard photons. Finally, we stress that despite 
the opposing trends in the temperature profiles shown in 
Figure for a fixed value of ^ the reflected spectrum is 
always softer for smaller Ecut- 

To produce the complete relativistic model relxill, 
we convolve spectra computed along the disk using 
xillver with the general-relativi stic ray tracing code 
relline (Dauser et al.|2010 2013). The model relxill 
self-consistently calculates the angle-dependent, rela¬ 
tivistic ionized reflection spectr a ([Garcia et al.|2014 ) that 
are used in fitting data. Figure [^shows reflected spectra 
(and the incident power-law spectra) for log ^ = 1 calcu¬ 
lated with relxill. For clarity, we show only two mod¬ 
els, for i5cut=300 keV and Acut=l MeV. These are the 
same as the corresponding two models shown in the top- 
left panel of Figure except that here the spectral fea¬ 
tures are blurred by Doppler and gravitational-redshift 
effects. 

The lower panel in Figure shows the ratio of the 
Ecut=300 keV model to the £'cut=l MeV model and also 
the ratio of the incident power-law spectra. In the NuS¬ 
TAR band, while the two models can be distinguished 
by considering either the power-law or the reflected com¬ 
ponent, it is the stronger curvature of the reflected com¬ 
ponent that is more telling. At lower energies, in the 
Suzaku band, the power-law is incapable of distinguish¬ 
ing between the two models because the continuum is 
essentially identical (apart from the normalization). On 
the other hand, the reflected component shows signifi¬ 
cant changes in spectral features, which sharply distin¬ 
guish the models. Thus, whether one considers NuSTAR 
alone, or NuSTAR plus Suzaku, it is possible to constrain 
Ecut at energies far above the bandpass of the detector, 
but in either case it is the reflected component that car¬ 
ries the greater weight. 

3. FITS TO SIMULATED NuSTAR DATA 

We first demonstrate that NuSTAR data alone can pro¬ 
vide strong constraints on Ecut- To this end, we sim¬ 
ulate data using appropriate response and background 
file^ assuming a circular extraction region with 60" ra¬ 
dius centered V off axis. We combine the data from 
both NuSTAR telescopes (focal plane modules FPMA 
and FPMB) and simulate spectra for an exposure time 
of 100 ks. 

We consider sources, both extragalactic and Galactic, 
with intensities/fluxes that cover the range nor mally en¬ 
countered in p ractice. Using the ISIS package (Houck & 
Denicola 2000), we grouped the data so that the mim- 
mum sighai-to-noise ratio for each bin was 5 to ensure the 
soundness of our analysis. We additionally grouped 
the data to ensure appropriate sampling, given that the 
resolution of the NuSTAR detectors is 0.4/0.9 keV at 
6/60 keV. Specifically, we further binned the simulated 

® http://www.nustar.caltech.edu/page/response_files 
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Figure 1. The left panels show xillver reflection spectra for log^ = 1 (top) and \og^ = 3 (bottom) and for nine values of Ecut (see 
legend in lower-right panel). For all models, the illu minating spectrum is a power law with index F = 2 and the iron abundance is solar. 
For simplicity we only show angle-averaged spectra < |Garcfa et al.|2014[ |. The right panels show the temperature profile in the slab. 
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Figure 2. Top: Relativistically blurred reflection spectra (solid 
lines) computed using relxill for log ^ = 1 and the respective 
power-law spectra (dashed lines) of the illuminating source with 
slope F = 2. The red curves are for F’cut=300 keV and the blue 
curves for F’cut = l MeV. Bottom: Ratio of the £^cut=300 keV spec¬ 
tra to the FJcut=l MeV spectra. The ratio of the reflection/power- 
law spectra is plotted as a solid/dashed line. The shaded box spans 
the NuSTAR bandpass (3—79 keV). 


below and above 20 keV, respectiveljjj 
In the simulations, we used the model Tbabs*relxill 
(Sect ion where Tbab s models the interstellar absorp¬ 
tion ( Wilihs et al. 2000). We simulated spectra for eleven 
values ot the high-energy cutoff ranging from 50 keV 
to 1 MeV. The other model parameters, which are the 
same for all simulations, are: the hydrogen column den¬ 


sity Nh = 5 X 10^^ cm“^; dimensionless spin parameter 
a = cJ/GM^ = 0.998; emissivity index q = 3; inclination 
of the disk z = 45 deg; inner disk radius i?in = i?iscoj 
power-law T = 2; ionization parameter log^ = 1; iron 
abundance Ape = 1; and the reflection fraction Rf = 3, 
which is defined as the ratio of the reflected to the power- 
law flux in the 20-40 keV band. Our calculations in this 
instance do not link the reflection fraction to the emissiv¬ 
ity index, which would require specifying the geometry 
of the reflector and the illuminating source. 

We emphasize that this Letter deals exclusively with 
relativistically-blurred reflection by ionized gas in the 
very inner portion of the accretion disk for which the 
reflection fraction can vary widely, with values ranging 
from a few tenths to > 20 for rapidly spinning black holes 
(Dauser et al. 2014). Reflection in this regime is not to 
be contused with reflection by distant, neutral material, 
which is observed in the spectra of many sour ces and is 
usually modeled using the pexrav model (e.g., Zdziarskij 
eral][l9^[200^ . 


We generated the sim ulated spectra using the fake it 
task in XSPEC v-12.8.2b ( Arnaud||1996 ), and we then fit- 


data, which have uniform bin widths of 0.04 keV over the 
entire energy band, to 0.2 keV and 0.4 keV for energies 


® We found that binning the data to a few samples per resolution 
element is important: Without this binning, the fitted values of 
Ecut were consistently biased below the input. This result suggests 
more generally that fitting oversampled data is not good practice. 
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Figure 3. Results of fitting simulated NuSTAR data. Top: Shown plotted are the values of Ecut returned by the fits compared to the 
values of Ecut used in producing the simulated data for three assumed values of source intensity (see legend). Fits using Model 1 (left 
panel) demonstrate that -Ecut can be accurately recovered for energies far above the 79 keV limit of the NuSTAR bandpass. Meantime, 
the very poor performance of Model 2 (right panel) is obvious. Error bars are at the 90% level of confidence. Bottom: Goodness of fit; 
note the poor performance of Model 2 at low to moderate values of Ecut for the brightest source. 


ted these spectra using two different models, with Ecut, 
r, Ape, normalization, and Rf free to vary: 

Model 1: Tbabs*relxill. We first employ the same 
model used to produce the simulated data in order to 
evaluate how well the parameters, -Ecut in particular, can 
be recovered. 

Model 2: Tbabs*highecut*relxill. In this case, 
the parameter Ecut in relxill is fixed at 1 MeV, while 
it is allowed to vary freely in highecut, a simple expo¬ 
nential convolution model acting on relxill. Model 2 
aims to test how accurately Ecut is recovered when the 
spectrum is described by modifying the high-energy re¬ 
gion of an otherwise incorrect reflection model (i.e., when 
ignoring the effects on the disk’s ionization structure de¬ 
scribed in §J^. 

Results ofmtting the simulated data are shown in Fig¬ 
ure The top panels show the fitted values of Ecut ver¬ 
sus the values used in simulating the data, with results 
for Model 1/2 in the left/right panel. The shaded regions 
indicate uncertainties at the 90% confidence level. The 
colors are keyed to source intensity. The quality of the 
fits is indicated in the lower panels. 

Focusing on the results for Model 1 (top-left panel) 
and the 10 mCrab source, which delivers 1.8 x 10® counts 
in the NuSTAR band, the recovered values of Ecut be¬ 
low 400 keV are of marginal quality, while above that 
energy they only provide lower limits on Ecut- How¬ 
ever, for source intensities of ^ 100 — 500 mCrab 
10^ — 10® counts), significantly better performance is 
achieved. We conclude that for high-quality data one can 
determine the high-energy cutoff with reasonable preci¬ 


sion using NuSTAR data alone. 

Turning to the results for Model 2 (top-right panel), 
one immediately sees the much poorer performance 
achieved in attempting to recover Ecut by simply fitting 
the curvature in the power-law continuum component us¬ 
ing the highecut model. As mentioned in Section!^ this 
was the standard approach used before Ecut was induded 
as a fit parameter in reflection models. The highecut 
model provides constraints of reasonable quality only for 
Ecut< 300 keV. Furthermore, for the brightest source 
considered (500 mCrab), Model 2 yields unsatisfactory 
results for -Ecut> 300 keV because not only is Ecut sig¬ 
nificantly over-predicted, the error bars are also mislead¬ 
ingly small. Interestingly, at lower energies the qualiW 
of the fits are worse for the brightest source (Figure 
bottom-right panel). Other model parameters are also 
poorly constrained, e.g., the Fe abundance is grossly un¬ 
derestimated, while both the power-law index and the 
reflection fraction Rf are over-predicted. 

Importantly, we note that Models 1 and 2 give the same 
quality of fit (with the exception of the 500 mCrab sim¬ 
ulation for low Ecut) so that a y^-test does not generally 
promote one model over the other. However, as we have 
demonstrated, Model 2 returns incorrect estimates of the 
parameters in all cases. In short, good statistical preci¬ 
sion does not ensure that the parameter values returned 
by a fit to a model are reliable or physically meaningful. 

4. INCLUDING LOW-ENERGY DATA 

Figure indicates that the low-energy component 
(E < 3 keV) of the reflection spectrum is more sensitive 
to the value of Ecut than the component in the NuSTAR 
































Constraining the high-energy cutoff 


5 



1000 

-M 

800 

Li- 

1 


> 

600 

a» 




■t-i 

400 

u 


LU 

200 


0 


1 mCrab 
10 mCrab 


0 


200 400 600 800 

Ecut (keV) - Real 


1000 


Figure 4. Including low-energy Suzaku data, Model 1 delivers extraordinary performance in constraining i?cut out to extreme energies 
compared to using NuSTAR data alone. The results shown here of fitting simulated data for a 10 mCrab source are superior to those 
obtained for a 500 mCrab source using only NuSTAR data (see the upper-left panel in Figure [^. 


band. Therefore, we now consider spectra comprised of 
Suzaku plus NuSTAR data in the expectation that this 
will significantly improve the constraints on Ecut, and 
the other parameters as well. The simulated low-energy 
data are prepared using Suzaku response files, the same 
set of model parameters, and the same source fluxes as 
those used in preparing the NuSTAR data. 

In the following, we consider only Model 1. Fitting 
the combined set of Suzaku (0.5-8 keV) and NuSTAR 
(3-79 keV) simulated spectra, we obtain the outstanding 
result shown in Figure]^ The fit for the 10 mCrab source 
is incomparably better than for the NuSTAR data alone 
(compare upper-left panel of Figure]^. The high-energy 
cutoff is strongly constrained for all values of Ecut con¬ 
sidered; even at 1 MeV, the error is only about 7%. This 
excellent performance of our reflection model is a conse¬ 
quence of the sensitivity of the reflection component to 
changes in the high-energy portion of the illuminating 
spectrum (Figure^. These simulations clearly demon¬ 
strate that fits to features in the relatively low-energy 
portion of the spectrum, where the instruments are most 
sensitive, can deliver accurate measurements of Ecut at 
energies far beyond the limit of the detector bandpass. 

Also shown in Figure Hare the results for a 1 mCrab 
source (2 x 10® counts). Remarkably, with the inclusion 
of low-energy data even ~ 10® counts are sufficient to 
recover Ecut with reasonable precision; the error bars for 
Ecut> 600 keV are ^ 30%. 


Finally, with two additional simulations we demon¬ 
strate that the reflection component is much more im¬ 
portant in constraining Ecut than the power law. For the 
same 1 mCrab source, and including Suzaku data. Fig¬ 
ure ^ compares the performance of a cutoff power-law 
model alone (cutoffpl) with a pure reflection model. 
These two limiting cases can be compared directly to the 
intermediate case of our adopted composite model (i.e., 
power law plus reflection) with Rf = 5 shown in gray in 
Figure The three cases considered together show that 
stronger reflection translates into a better constraint on 
A^cut- 

5. DISCUSSION AND CONCLUSIONS 

We have demonstrated that our reflection model 
relxill can constrain the physically important coronal 
parameter Ecut for spectra with intensities like those en¬ 
countered in observations of AGN. We have furthermore 
shown that attempting to constrain Ecut by imposing a 
cutoff that is external to the reflection model yields rela¬ 
tively poor constraints on Ecut- Such a simplified model 
inaccurately represents the observed spectrum and in¬ 
troduces strong biases in the values inferred for the key 
physical parameters. 

We analyzed two simulated data sets using relxill; 
(1) With NuSTAR data alone (3-79 keV), Ecut> 400 keV 
can only be constrained if the spectrum contains > 
10^ counts, which is infeasible for most observations of 
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Figure 5. Results simulated for a 1 mCrab source demonstrating that a pure reflection model (i.e., no power-law continuum) shown in 
red is far superior to an e-folded power-law component (i.e., no reflection) in constraining -Ecut* 


AGN. (2) However, our simulations indicate that for a 
simultaneous fit to Suzaku and NuSTAR data (0.5-79 
keV), one can constrain Acut for a 1 mCrab source at 
energies far above the limit of the detector bandpass. 

Our results lead to two basic conclusions: 

1. In constraining the cutoff energy, misleading re¬ 
sults are obtained if the power-law component is 
cut off externally. It is therefore essential to em¬ 
ploy a detailed and accurate model of the reflection 
spectrum that includes Ecut as an internal fit pa¬ 
rameter. 

2. Because the cutoff energy most strongly affects the 
low-energy part of the reflection spectrum, to ob¬ 
tain useful constraints on Ecut for AGN it is crucial 
to combine NuSTAR data with low-energy (e.g., 
Suzaku or XMM-Newton) data. 

This limited study provides only a very cursory explo¬ 
ration of how proper modeling of a reflection spectrum 
yields useful physical constraints on the high-energy cut¬ 
off out to extreme energies. A thorough exploration of 
the subject is well beyond the scope of this Letter. For 
example, we illustrate the effects on the reflected spec¬ 
trum of varying Ecut for the cases of log^ = 1 and 3 
(Figure while only presenting simulations for the for¬ 
mer case. However, one expects that the effects will be 
quite different for other choices of the ionization param¬ 
eter. Moreover, the values of other parameters, such as 


the reflection fraction, can affect the statistical signifi¬ 
cance of the value returned for Ecut- 
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